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ABSTRACT 

The origin of the extragalactic gamma-ray background is a pressing cosmological mystery. The Fermi 
Gamma-Ray Space Telescope has recently measured the intensity and spectrum of this background; 
both are substantially different from previous measurements. We present a novel calculation of the 
gamma-ray background from normal star-forming galaxies. Contrary to longstanding expectations, we 
find that numerous but individually faint normal galaxies may comprise the bulk of the Fermi signal, 
rather than rare but intrinsically bright active galaxies. This result has wide-ranging implications, 
including: the possibility to probe the cosmic star-formation history with gamma rays; the ability to 
infer the cosmological evolution of cosmic rays and galactic magnetic fields; and an increased likelihood 
to identify subdominant components from rare sources (e.g., dark matter clumps) through their large 
anisotropy. 

Subject headings: gamma rays: diffuse background — gamma rays: galaxies — cosmic rays — galaxies: 
star formation 



1. INTRODUCTION 

The Fermi Gamma-Ray Space Telescope has unveiled 
the high-energy cosmos with unprecedented clarity and 
depth. The gamma-ray sky has been known (e.g., 
iHunter et all 119971: iSreekumar et al.lll998| ) to be domi- 
nated by diffuse emission from the Galactic plane, while 
at high Galactic latitudes a diffuse extragalactic gamma- 
ray background (EGB) has an important, and at some 
energies dominant, contribution. However, before Fermi, 
the processes dominating the diffuse emission from the 
Galaxy, especially above 1 GeV, were unclear and highly 
debated-cf. discussion on the Ge V excess reported by the 
Ener getic Gamma-ray E xperiment Telescope (EGRET) 
(e.g.. IStrong et all 120001 and re ferences therein). Fermi 
has clarified (fAbd o et al. 2009a) that the dominant emis- 
sion mechanism is cosmic-ray interactions with interstel- 
lar gas, which leads to gamma rays mostly from pion 
decay in fligh t, i.e., p CI + Pi sm — > ppn° then ir° — > 77 
(jSteckedll97in . 

Moreover, the Fermi EGB differs from previous 
EGRET estimates: the intensity is fainter and the spec- 
trum steeper, consistent with a power law of spectral in- 
dex 2.41 ± 0.05 and integrated inten sity J(> lOOMeV) = 
(1.03 ±0.17) x lO^cnr^s^sr" 1 (|Abdo et al.l I2010bf) . 
Our theoretical understanding of the EGB must therefore 
be substantially revised in light of the new and smaller 
Fermi signal. 

Pioneering studies investigating the origin of the 
EGB first considered the collective emission form star- 
forming galaxies (like the Milky Way), but found 

bdfields@illinois.edu 
1 Einstein Fellow 



this to give a sma ll EGB signal (IStrong et all 119761 : 
IBignami et aT1ll979l ). IPavlidou fc Fields! (<2002f )first in- 
corporated observations of the c osmic star- formation 
rate, while iProdanovic fc Fields! (|2006l ) made the first 
estimates of the pionic contribution from star-forming 
galaxies. In both cases, the predicted intensity was 
below the then- measured EGB. Blazars, the bright- 
est extr agalactic sources, have be en the favored can- 
didates (jStecker fc Salamonl 119961 ). However, subse- 
quent esJimates__of their contribution to the EGB 
(e.g.. IMiicke fc Pohll 120001 : IChiang fc Mukheried Il998t 
iNarumoto fc TotanTl2007t iDermerl l2007f) have consis- 
tently fallen short. Fermi point-source observations sug- 
gest that unresolved blazars contribute at most ~ 23% 
of the EGB : and t hus the mystery has become acute 
(jAbdo et allboiOdn . 

Here we present a more realistic model for the EGB 
from star-forming galaxies, constructed to use as much as 
possible of our substantially improved multiwavelength 
observational understanding of these sources, and isolat- 
ing the signal from normal star formation (as opposed to 
starburst galaxies). 

2. FORMALISM 

The EGB intensity is an integral of the gamma-ray 
luminosity density C~ t (emissivity) over the line-of-sight 
to the cosmic horizon 
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with H(z) the Hubble function. We 
adopt a ACDM cosmology, with a Hubble parameter 
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Hq = 71 km s _1 Mpc -1 , and cosmological constant and 
matter density parameters Qm = 0.3 and £1a = 0.7 re- 
spectively. 

The cosmological inputs to eq. ([T]) are well-determined, 
so the EGB entirely hinges on the luminosity density £ 7 
and its connection to cosmic star formation. In this study 
we construct, for the first time, a gamma-ray luminosity 
function (distribution of sources by gamma-ray luminos- 
ity as a function of rcdshift) for normal galaxies. The 
luminosity function is used to obtain the EGB intensity, 
by integrating over luminosities and redshift along pho- 
ton paths. 

The star-forming luminosity density due to pionic 
emission follows from the gamma-ray luminosity per star- 
forming galaxy L~. 



(2) 



where n 7 is the comoving number density of gamma-ray- 
luminous galaxies. The average is taken over the dis- 
tribution of star- forming galaxy properties at redshift z, 
i.e., the gamma-ray luminosity function which gives the 
comoving number n 7 (L 7 ) c?L 7 of star- forming galaxies 
with luminosity in the range (L 7 ,i 7 + c?L 7 ); out to a 
maximum L max determined below. 

A galaxy's pionic gamma-ray flux scales with the 
cosmic-ray flux (projectiles) and the amount of interstel- 
lar gas (targets) in the galaxy. Specifically, a galaxy's 
rest-frame pionic gamma-ray luminosity (photon counts 
per unit time) is given by 
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with E em the photon energy in the emitting frame. Here 
the pionic gamma-ray production rate per interstellar H- 
atom is r„.o_ ! . 77 (£ , cm ) = (<£> p dovo_> 77 /<iE cm ) oc $ p and 
is proportional to t he galaxy's volume-averaged cosmic - 
ray proton flux $ p HPonHll994t IPersic k, Rephaelill20T(l ; 
the cross-section da^o^ rll /dE cm is understood to include 
effects of pion multiplicity and of helium and heavier el- 
ements in the cosmic rays and interstellar medium. The 
factor Af p — J nn dVi sm gives the number of hydrogen 
atoms in the galaxy's interstellar medium, summed over 
all states-molecular, atomic, and ionized. This term 
is fixed by the galaxy's total interstellar gas mass via 
TVh = XuMgas/nip, where Ah w 0.70 is the hydrogen 
mass fraction. 

Since cosmic rays are thought to be predominantly 
accelerated in sites associated with massive star for- 
mation (supernova remnants, massiv e stellar winds, 
pulsars, see e.g., ILacki et all ()2010D and references 
therein), their flux in eq. © should scale with the 
star formation rate (SFR) in the galaxy. Observa- 
tions by Fermi and TeV telescopes HESS and VER- 
ITAS have confirmed this expec tation through ob ser- 
vations of star- fo r ming galaxies (lAbdo et al.l I2010cllallel ; 
lAcero et al.|[2009t IVERITAS Collaboration et all 120091 ) . 
for which both SFR and total gas mass were indepen- 
dently known. Thus the cosmic-ray flux is set by the 
competition of particle acceleration and losses that we 
assume are dominated by escape, as they are in the 



Milky Way. These self-regulating processes tend toward 
an equilibrium & p oc A csc ^, where ip is the star forma- 
tion rate and A esc is the escape pathlength (assumed con- 
stant). We thus adopt the scaling 

L 1 cx M gas ip (5) 

of gamma-ray luminosi ty with a galaxy's gas mass and 
star- formation rate ip (iPavlidou fc F ields 200l|) . 

The scaling in eq. ([5]) represents a central ansatz re- 
garding gamma-ray production in Milky- Way-like galax- 
ies whose cosmic-ray losses are dominated by escape. 
This physically-motivated relation now has support from 
Fermi observation of star-forming galaxies for which the 
product M gas ip ranges over several orders of magnitude. 
A second ansatz in our model is our adoption of a univer- 
sal cosmic-ray spectrum of index s cr = 2.75. At energies 
away from the peak the gamma-ray spectrum has the 
same index: s 7 .i ia d = s CI . Even though both approxima- 
tions are certainly simplifying, they represent important 
benchmark cases, against which more sophisticated mod- 
els can be tested. 

The scaling law of eq. (O allows us to determine the 
gamma-ray output of any star-forming galaxy, but only 
once we normalize our results to a system in which the 
cosmic-ray and/or gamma-ray properties are known. We 
choose to normalize to the Milky Way, where the local 
cosmic-ray flux is well-measured, and the global star- 
formation rate is also known. We assume that the ra- 
tio of cosmic-ray flux to star-formation rate should be 
constant for all normal galaxies, i.e., that 



7T U — f 77 



rMW 

1 -rrO. 



77 



^cr.MW -RSN.MW V'MW 



(G) 



This scaling encodes the longstanding notion that super- 
nova remnants accelera te hadronic cosmic rays (leading 
to pionic emission, e.g. JAbdo et al.ll2009cf) . 
We thus find that for normal galaxies, 

^(M gllSl lp)—X H L w , 
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We adopt the photon emission per hydrogen atom 
derived from Fermi diffuse Galactic observations at 
medium latitudes for photons > 100 MeV: T^™ 77 = 

2.0 x 10~ 25 s^H-atonT 1 (jAbdo et al.l l2009b[) . This 
value represents a large-scale (~kpc) spatial averaging 
of Galactic cosmic-ray properties, which is appropriate 
for our global calculation. 

We can now write the gamma-ray luminosity density 
for normal galaxies 



£y(E C ni, z) — Ah 



r^ 77 (£ cm ) (M gas (z)) . 



P*(z) (8) 



in terms of the cosmic star-formation rate pi,(z) — 
(^rigai) . We define a mean interstellar gas mass as 

(Mgas (*)> EE i%^> = (*W%al) (g) 

The assumption that losses are escape-dominated and 
uniform across galaxies can only be approximately valid 
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at best. This is a major uncertainty in our model, which 
will benefit from future data on the EGB and resolved 
galaxies. For example, even the leaky-box model can 
generalize eq. ([5]) to L 7 oc A csc ^M sas ; variations in the 
energy dependence of the escape length A osc would also 
change the cosmic-ray and photon spectral indices which 
we take as universal. 

Indeed, Fermi observations of the Large Magellanic 
Cloud suggest that cosm ic-ray confinement and prop- 
agation are non-trivial (|Abdo et al.l l2010d ). More- 
over, starburst galaxies show very high cosmic-ray in- 
tensities within small volumes where inelastic colli- 
sions compete with, and someti mes dominate, out- 
flows to regulate cosmic-ray losses dPaglione et al .11X9961; 

T nr.\A of ol I lOm nl IT^vt-qo If nil 1000/11 Tln^m^o^n of ol I 



Lacki et all 120101 iTorres et all 12004 [Thompson et al.l 
20071 iPersic fc Rephaeli I2010t fstecken l2007f ) . For this 
reason, eq. ((7} provides a rough description of normal 
escape-dominated galaxies only; we do not expect it to 
hold for starburst galaxies, which we will exclude below. 

2.1. Gas Mass and Star- Formation Rate 

We can infer a galaxy's interstellar gas mass at a given 
star-formatio n rate, via the well-established Kennicutt- 
Schmidt law (|Schmidtlll95l lKennicuttJll998t ). The sur- 
face densities for star formation and gas are found to 
be correlated via i^/Moyr^kpc" 2 = (2.5 ± 0.7) x 
lO^Egas/Mopc- 2 )* with x = 1.4 ± 0.15. Both 
normal and starburst galaxies follow this correlation, 
but normal galaxies populate £*, no rm a l < £™normai = 
0.4 M© yr -1 kpc~ 2 while starbursts occupy the opposite 
regime (|Kennicuttill998f ). 

To recover a relationship between the global galac- 
tic star-formation rate tp = 7rrj isk S Jr and gas mass 



7rr 



disk 



_£ gas , requires a galactic disks scale length 
fVtisk- We take rnisk = 18.9 kpc/(l + z); observations 
(|Erb e t al. 2006) indicate that this choice is uncorrelated 
with galaxy mass and star-formation rate at a fixed z. 
Combining this with the Kennicutt-Schmidt law, we find 



M gas (V, z) = 2.8 x 1O 9 M (1 + z) 



4> 



1 M yr-i 



(10) 

where /3 = 2(1 - l/x) = 0.571 and w = l/x = 0.714. In 
our model, normal galaxies extend to 

V> < VC a r X ma iM = 7rr disk (z) 2 I]™* r mal ~ _^L_M /vr 

(11) 

this cutoff becomes important for z > 1. 

Combining eqs. (fi"0)) and (0, we can express a galaxy's 
gamma-ray luminosity in terms of its SFR and redshift: 



L 7 fcM« (l + z)- f3 
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Available Fermi data on resolved z ks star-forming 
galaxi es are consistent with this scaling ()Abdo et al.l 



2.2. Gamma-Ray Luminosity Function 

A galaxy's luminosity in the Ha line provides a well- 
established tracer of star formation rate: Lhq/1-26 x 



10 34 W = ip/l Mq yr" 1 (|Hopkinsll200l . The Ha distri- 
bution of galaxies (luminosity function) is related to the 
cosmic star-formation rate density via p*(z) — {ipn) oc 
/ I/Ha Ti(L-a a ,z) dLn a , which by now is well-measured 
as a function of redshift both by Ha and by other meth- 
ods, so that its shape is well-determined observationally. 

The gamma-ray luminosity density is a different mo- 
ment of the Ha luminosity function C 1 {z) — (ipn) oc 

(1 + z) - ' 3 J Llia £jja n (Lna,z) dL-Ra via our scaling laws 
above. We see that the gamma-ray luminosity distribu- 
tion (i.e., luminosity function) at a given redshift follows 
directly from the distribution of star-formation rates, 
as traced by Ha. Since we consider here only normal 
galaxies, we include only galaxies with ip < Vernal 
(eq. [TT]) which sets corresponding limits L^ia < ^HoS 
and L 7 < in eq. ©. 

Current data on the Ha luminosity function can be 
fit to a Schechter function, of the form n(L, z) dL = 
(L/L*)~ a cxp (— L/L*) dL/L* where L — L^ a . 
Present data are consistent with the value a = 1.43, 
found for z = 0, persisting for all redshifts. Data also fix 
the z = values ra»(0) = 1. x 10~ 3 Mpc" 3 and £ *(0) = 
9.5 x 10 34 W, for h = 0.71 (jNakamura et al.ll200l . This 
corresponds to a star-formation rate ip(L r ) — 7.5 M /yr. 

However, observations currently do not give unambigu- 
ous solutions for the other two parameters, the char- 
acteristic comoving density of star-form ing galaxies n * 
and the characteristic Ha luminosity L r (Hopkins 2004) . 
Moreover, it is unclear whether and how each parameter 
evolves with redshift. 

Two limiting cases bracket the possible behaviors of the 
Ha luminosity function and thus of cosmic star forma- 
tion. In the case of pure luminosity evolution the comov- 
ing density of stars is fixed, n* = const independent of z, 
and all redshift evolution lies in I* = L sr {z). Conversely, 
pure density evolution places the redshift evolution in 
the density scale n*(z) while setting L» = const. In each 
of these two limits, we can find the redshift dependence 
of the parameters via the requirement (Ln) oc p*. The 
redshift history of co smic-star formation is well-known 
(jHoriuchi et alJl2009f) . and we encode this in the dimen- 
sionless "shape" function 



S(z) = p4z)/p4Q). 



(13) 



We then have L*(x)/I*(0) = S(z) in the case of pure 
luminosity evolution, and n.*(z)/n*(0) = S(z) in the case 
of pure density evolution. 

At a given redshift, eq. © gives the scaling (M gas ) = 
(M gas ipn) /p± oc (M gas tpn) / S (z) . For our Ha luminos- 
ity function and Kennicutt-Schmidt relation we find a 
local value of (M gas )2 =0 = 6.8 x 10 9 M . At other red- 
shifts we have (M gas ipn) oc (1 + z)^ 13 L^(z) LU S(z), and so 
(M gas ) oc (1 + z)- L*(z) u . Thus for the pure luminos- 
ity case -L* oc S(z), we find that the gas mass strongly 
evolves as (M gas ) oc (1 + z)~° S(z) u , in response to the 
strongly changing SFR. Consequently, the factor of 10 
rise in cosmic star-formation at z ~ 1 implies a net 
gamma-ray luminosity increase of a factor ~ 30. On the 
other hand, in the pure density evolution case, galaxy 
SFRs are constant, L*(z) = const, so that the mean gas 
mass (Mg as ) oc (1 + z)~P actually decreases with red- 
shift, while the comoving number of star-forming galax- 
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ies increases, but the net enhancement at high redshift is 
smaller than in the pure luminosity evolution case. This 
key difference leads to the factor ~ 4 between the EGB 
predictions for the pure density and pure luminosity evo- 
lution cases seen in Figure 1. 

3. RESULTS AND CONCLUSIONS 

Our full numerical calculation uses a Milky Way 
pionic source spect r um w hose shape is derived from 
Pfrommcr & Enfilin (2004), calibrated to observations 
by normalizing the > 100 MeV photon emission per hy- 
drogen atom to the Fermi result at intermediate Galac- 
tic lat itudes (lAbdo et all l2009bl) . The cosmic SFR is 
from iHoriuchi et all (j2009fl . ~For the Milky Way SFR, 
used to normalize the co smic-ray flux/SFR rat i o, we 
use the recent estimate of [Robitaille fc Whitnevl ()2010l ) 
(V'mw = 1-^0 /y r j a factor of 3 lower than earlier work). 

Figure 1 shows our results for the normal galaxy 
contribution to the EGB. We plot predictions for the 
limiting cases of pure luminosity and of pure den- 
sity evolution. The uncertainties in the model inputs, 
summed in quadrature, propagate into the displayed er- 
ror band that applies to each curve, which we estimate 
to be a factor of 10 ±03 , res ulting from uncert ainties 
of: 30% in pionic emissivity (|Abdo et all l2009bj) . 40% 
in the normalization of the G alactic star-formation rate 
(|Robitaille fc Whitnevl I2010D. 40% in the cosmic star- 
formation rates (|Horiuchi et al.l 120091 ) . and 25% in the 
luminosity scaling in eq. (j9|). The true systematic uncer- 
tainty would also reflect the idealizations in our model 
(universal cosmic-ray spectra and confinement). These 
errors are hard to estimate but in any case imply that 
the uncertainty range in Figure 1 is a lower bound to the 
error budget. 

Within errors, our predictions for both limiting mod- 
els fall at or below the level of the Fermi data, where 
the data seem to support the pure luminosity evolu- 
tion case that explains nearly the entire signal. Com- 
paring central values, this model gives ~ 50% of the 
Fermi EGB < 10 GeV. Thus, unresolved normal galax- 
ies make a substantial and likely dominant contribution 
to the observed EGB, without overpredicting the signal. 
Even the pure density evolution case accounts for a min- 
imum of 20% of the EGB around 0.3 GeV; this provides 
a lower limit to the normal-galaxy signal Thus, any 
other EGB sources dStecker fc Salamonl Il996t iDermerl 
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120071 : lAbdo et al.M2010dt iFaucher-Giguere fc Loebll2010j ) 
must contribute no more than the remaining 80% of the 
data. Indeed, the LAT team upper limit to the blazar 
EGB contributi on shown in Figur e 1 is comparable to 
our lower limit (|Abdo et al.ll2010dH 

The spectral shapes of the two limiting cases are very 
similar: the peak in E 2 dI/dE lies at ~ 0.3 GeV because 
the bulk of the signal comes from z ~ 1. These mod- 
els predict that the EGB turns over for E < 0.3 GeV, a 
testable prediction of our model. For hadronic emission, 
the high-energy spectral index is the same as the under- 
lying proton spectral index, here s 7i had = Scr = 2.75; this 
is somewhat steeper than the Fermi single-power-law fit 
s bs = 2.41 ±0.05. Consequently, our predictions at high 
energies (> 10 GeV) fall below the data. If normal galax- 
ies had a distribution of cosmic-ray spectral indices, the 
resulting EGB spectrum would steepen at high energies 
where the hardest sources would dominate, developing a 
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Fig. 1. — The normal galaxy contribution to the extragalac- 
tic gamma-ray background. Curves represent two limiting cases 
of cosmic star formation: pure luminosity and pure density evolu- 
tion. Colored error band illustrates the factor 10 ±0 ' 3 uncertainty 
in the norm alization of both theory curves. Fermi data are from 
lAbdo et al. 2010 E|), Do tted line: unresolved blazar EGB upper 
limit IIAbdo et al.ll2010dl) . 

convex tail. Indeed, the Fermi EGB data suggests a slight 
flattening of slope around E > 10 GeV, which might hint 
at such a transition. 

A galaxy with characteristic Ha luminosity has 
L*{> lOOMeV) = 1.4 x 10 43 s" 1 . Such objects have 

flux F if they lie at distances r* = (L*/47rF) 1/2 

I 1 Mpc (10~ 9 cm- 2 s^/F) 1 / 2 . Thus Fermi should 
eventually resolve 



N(> F) ~4nrl n.(0)/3 = 



io- 



3/2 



(14) 

normal galaxies, consistent with 2-3 de tections to 
date ( th e LMC, SMC, and perhaps M31: lAbdo et all 
I2010ct iThe Fermi/LAT collaboration fc Abdol 120101: 
lOgelman et~alll2010D . 

Our results do not account for starburst galaxies, 
nor for inverse- Compton emission from any star-forming 
galaxies; these must contribute to the star-forming EGB, 
and could have hard spectra dominating > 10 GeV. We 
have also neglected gamma-ray attenuation by extra- 
galactic background light (important at E > 30 GeV; 
e.g.. lStecker et aill2006l and references therein). We will 
address these issues in future work. 

The amplitude and configuration of magnetic fields in a 
galaxy have an additional effect on the scaling of cosmic- 
ray flux with SFR. Confirmation that normal galaxies 
comprise the bulk of the Fermi signal would constitute 
a unique probe of the evolution of these magnetic fields 
between the redshift of peak star formation and today. 

Because of their ubiquity, normal galaxies produce the 
smallest anisotropics in the EGB, far less than blazars or 
other proposed sources. Thus, by studying the observed 
EGB anisotropy as a function of energy, it may be pos- 
sible to disentangle the spectrum and amplitude of the 
normal-galaxy contribut ion from that of other sources 
(|Ando fc Pavlidoul 12009ft . Moreover, because normal 
galaxies seem to dominate the Fermi EGB, their small 
contribution to anisotropies will fortuitously provide the 
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optimal chance of findi ng smaller and more exotic source s 
in the observed signal (Siegal-Gaskins & Pavlidou 2009). 
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